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Abstract 

A  single  reversible  solid  oxide  fuel  cell  (RSOFC)  can  accomplish  two  functions:  (1)  as  a  solid  oxide  steam  electrolyzer  (SOSE)  for  hydrogen 
production  and  (2)  as  a  solid  oxide  fuel  cell  (SOFC)  for  power  generation.  An  electrochemical  model  was  developed  to  study  the  performance  of  an 
RSOFC  based  on  a  proton-conducting  electrolyte  (RSOFC-H).  In  both  SOSE  and  SOFC  modes,  the  hydrogen  electrode- supported  configuration 
was  identified  as  the  most  favorable  design  to  achieve  high  energy  conversion  efficiency  of  RSOFC-H.  For  comparison,  in  a  previous  study  on 
conventional  RSOFC  based  on  an  oxygen  ion-conducting  electrolyte  (RSOFC-O),  the  hydrogen  electrode- supported  configuration  was  found  to  be 
favorable  in  the  SOFC  mode  but  such  configuration  would  cause  high  concentration  overpotential  in  the  SOSE  mode.  Thus,  the  oxygen  electrode- 
supported  configuration  was  desirable  for  RSOFC-O  operating  in  the  SOSE  mode.  The  results  obtained  in  this  study  show  that  RSOFC-H  has  a 
natural  advantage  over  RSOFC-O  in  terms  of  structural  design.  The  modeling  study  signifies  the  difference  between  RSOFC-H  and  RSOFC-O 
and  can  serve  as  a  useful  tool  for  further  design  optimization. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A  reversible  solid  oxide  fuel  cell  (RSOFC)  can  be  operated  in 
the  SOSE  mode  to  produce  hydrogen  fuel  and  the  same  RSOFC 
can  later  consume  the  hydrogen  fuel  to  generate  electricity  in 
the  SOFC  mode  [1-3].  Due  to  the  site-specific  and  intermit¬ 
tent  characteristics  of  renewable  resources,  such  as  solar  and 
wind,  hydrogen  is  a  potentially  useful  energy  carrier  to  sup¬ 
port  renewable  utilization.  The  capability  of  performing  dual 
function  using  a  single  cell  can  thus  make  RSOFC  potentially 
economically  sound  to  better  utilize  renewable  resources. 

Several  experimental  studies  on  conventional  RSOFC  based 
on  oxygen  ion-conducting  electrolyte  (RSOFC-O)  have  been 
performed  to  investigate  the  performance  characteristics  of 
RSOFC  in  both  working  modes  [4-6].  Recently,  a  modeling 
study  on  the  concentration  overpotentials  of  an  RSOFC-O  has 
been  conducted  [2] .  The  transport  and  reactions  involved  in  the 
two  working  modes  of  the  RSOFC-O  are  illustrated  in  Fig.  la 
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and  b,  respectively.  It  is  found  that  a  hydrogen  electrode  sup¬ 
port  is  favorable  to  the  SOFC  mode  while  the  oxygen  electrode 
support  is  favorable  to  the  SOSE  mode  [2].  The  difference  in 
the  optimal  designs  is  due  to  the  different  mass  transport  char¬ 
acteristics  of  gaseous  species  in  the  porous  electrodes  in  the 
SOSE  and  SOFC  modes.  This  study  indicates  that  RSOFC-O 
has  fundamental  problems.  No  matter  which  electrode  support 
(hydrogen  or  oxygen)  is  selected,  the  RSOFC-O  will  exhibit  high 
concentration  overpotential  in  one  of  the  two  working  modes. 

Alternatively,  an  RSOFC  can  be  built  with  a  proton¬ 
conducting  electrolyte  (RSOFC-H).  Although  some  experi¬ 
mental  studies  have  been  accomplished  to  characterize  the 
performance  of  RSOFC-H  in  the  SOSE  mode  and  SOFC  mode 
[7-10],  the  present  literature  lacks  detailed  modeling  analy¬ 
ses  on  RSOFC-H.  Modeling  studies  are  important  for  better 
understanding  of  the  working  mechanisms  as  well  as  for  design 
optimization.  Under  the  same  operating  conditions  (i.e.  temper¬ 
ature,  current  density),  an  RSOFC-H  working  in  either  SOFC 
or  SOSE  mode  has  identical  Ohmic  overpotentials.  In  addi¬ 
tion,  the  exchange  current  densities  for  SOSE  and  SOFC  are 
identical.  Therefore,  the  activation  overpotentials  of  SOFC  and 
SOSE  can  be  considered  identical  as  well  using  the  simplified 
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Fig.  1 .  Working  mechanisms  of  RSOFC-O:  (a)  SOSE  mode  and  (b)  SOFC  mode. 


Butler-volmer  equation  (omitting  the  reactant/product  con¬ 
centrations).  It  is  noted  that  the  activation  overpotential  is 
independent  of  the  electrode  thickness,  which  is  the  focus  of 
the  present  study.  These  overpotentials  can  be  calculated  with 
the  electrochemical  model  recently  developed  [11].  On  the  other 
hand,  the  concentration  overpotentials  are  different  between  the 
SOSE  mode  and  the  SOFC  mode  due  to  different  gas  transport 
characteristics  in  the  porous  electrodes.  Therefore,  the  concen¬ 
tration  overpotentials  are  the  sole  factor  responsible  for  different 
RSOFC-H  current-voltage  ( J-V)  characteristics  between  the 
SOSE  and  SOFC  modes.  For  this  reason,  this  electrochemi¬ 
cal  modeling  study  on  RSOFC-H  focuses  on  the  concentration 
overpotentials  in  both  SOSE  and  SOFC  modes. 

2.  Model  development 

2.7.  Operating  mechanisms 

The  fundamental  working  mechanisms  of  an  RSOFC-H 
working  in  the  SOSE  mode  are  illustrated  in  Fig.  2a.  Steam  is  fed 
to  the  porous  oxygen  electrode  and  transported  to  the  interface 
between  the  oxygen  electrode  and  the  electrolyte  where  electro¬ 
chemical  reactions  take  place.  The  important  electrochemical 
reaction  steps  include  (1)  adsorption  of  steam  onto  the  surface 
of  electrode  particles,  (2)  electrochemically  splitting  of  steam  to 
form  oxygen,  electrons  and  protons,  (3)  desorption  of  oxygen 
produced,  and  (4)  incorporation  of  proton  into  the  electrolyte. 


The  protons  are  transported  through  the  dense  electrolyte  to 
the  porous  hydrogen  electrode  where  they  are  combined  with 
electrons  to  form  H2.  It  should  be  mentioned  that  the  consump¬ 
tion  of  H2O  and  generation  of  O2  will  result  in  concentration 
gradients  of  these  two  species  in  the  porous  oxygen  electrode. 
In  addition,  the  molar  production  rate  of  O2  is  only  half  of 
the  molar  consumption  rate  of  H2O,  thus,  the  pressure  at  the 
electrode-electrolyte  interface  is  lower  than  that  at  the  surface 
of  the  electrode.  Therefore,  the  transport  of  H2O  and  O2  is  by 
means  of  both  diffusion  and  permeation.  Since  H2  is  the  only 
gas  in  the  porous  hydrogen  electrode,  it  is  transported  by  means 
of  permeation. 

The  working  mechanisms  of  RSOFC-H  in  the  SOFC  mode 
are  illustrated  in  Fig.  2b.  Similar  to  the  SOSE  mode,  H2  (as  a 
single  gas)  is  transported  in  the  porous  hydrogen  electrode  by 
means  of  permeation,  while  the  O2  and  H2O  are  transported 
in  the  porous  oxygen  electrode  by  means  of  both  diffusion  and 
permeation. 

2.2.  RSOFC-H  cell  potentials 

In  the  SOSE  mode,  the  J-V  characteristics  of  an  RSOFC-H 
can  be  expressed  as: 

ySOSE_£  ^actjn2+^actj  o2  +  ^Ohm  +  ^conc^  +  ^conc^  W 
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Fig.  2.  Working  mechanisms  of  RSOFC-H:  (a)  SOSE  mode  and  (b)  SOFC  mode. 
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where  E  is  the  equilibrium  potential;  ^act,H2  anc*  *7act,o2  stand 
for  the  activation  overpotentials  of  the  hydrogen  electrode  and 
oxygen  electrode,  respectively;  r]Ohm  represents  the  Ohmic 
overpotential;  and  0con^u2  anc^  1lcon^o2  are  ^e  concentration 
overpotentials  at  the  hydrogen  electrode  and  oxygen  electrode, 
respectively. 

In  the  SOFC  mode,  the  operating  potential  (yS0FC)  of  an 
RSOFC-H  can  be  determined  as: 


V' 


SOFC 


—  E  7?act,H2  ^act,02  ^Ohm 


SOFC 
I  conc,H2 


SOFC 
I  conc,02 


(2) 


where  Oconc^  anc*  ^concC02  stand  f°r  the  concentration  overpo¬ 
tentials  of  the  two  electrodes  working  in  the  SOFC  mode. 

The  equilibrium  potential  E  can  be  determined  with  the 
Nemst  equation  [12].  The  Ohmic  (?7ohm)  and  activation  over¬ 
potentials  (?7act,H2  and  ^act,o2)  as  well  as  the  concentration 
overpotentials  in  the  SOFC  mode  (0comk2  anc*  ^concEo2)  can 
be  determined  with  an  electrochemical  model  recently  devel¬ 
oped  [11].  The  concentration  overpotentials  in  the  SOSE  mode 
(^concCH2  an<^  ??concCo2)  can  determined  as  described  in  the 
subsequent  section. 


2.3.  Concentration  overpotentials 


In  this  section,  the  procedures  to  calculate  the  important  con¬ 
centration  overpotentials  of  RSOFC-H  in  the  SOSE  mode  are 
developed.  The  concentration  overpotentials  represent  the  resis¬ 
tance  of  the  porous  electrode  to  the  transport  of  gaseous  species 
[13].  According  to  their  definitions,  the  concentration  overpo¬ 
tentials  of  an  RSOFC-H  in  the  SOSE  mode  can  be  written  as: 


SOSE 

A/conc,H2 


and 


SOSE 

A/conc,02 


57L in  ( (/V/2<o) 

2F  {iP°of2Pko) 


(3) 


(4) 


where  P^2,  Pq2,  and  Ppj  Q  represent  pressure  of  hydrogen  at 
the  hydrogen  electrode-electrolyte  interface,  partial  pressure  of 
oxygen  at  the  oxygen  electrode-electrolyte  interface,  and  partial 
pressure  of  steam  at  the  oxygen  electrode-electrolyte  interface, 
respectively;  F  is  the  Faraday  constant;  R  is  the  ideal  gas  con¬ 
stant;  the  superscript  0  represents  the  surface  of  an  electrode. 

In  a  steady  state,  the  transport  of  each  participating  compo¬ 
nent  is  determined  by  the  local  conservation  of  mass  [14]. 


WNi  =  0 


(5) 


where  Ni  is  the  mass  transfer  of  species  i.  H2  is  the  only  com¬ 
ponent  present  in  the  porous  hydrogen  electrode.  The  transport 
of  H2  by  means  of  permeation  can  be  determined  with  Darcy’s 
law  [15,16]. 


Afii2  =  - 


PH2  Ba 
RT/x  h2 


VPh2 


(6) 


where  /xh2  is  the  dynamic  viscosity  of  H2  and  Bg  is  the  per¬ 
meation  coefficient,  which  can  be  determined  by  a  sixth-order 
polynomial  function  developed  by  Todd  and  Young  [17]  and  the 
Kozeny-Carman  relationship  [15],  respectively. 

Substituting  Eq.  (6)  into  Eq.  (5),  the  controlling  equation  for 
H2  permeation  in  the  porous  hydrogen  electrode  can  be  obtained. 


3  I"  Ph9  Bq  3Ph9 

—  2  g  — ^  =0  (7) 

dx  l RT/iu2  dx 

where  v  is  the  normal  distance  measured  from  the  electrode 
surface  and  is  applicable  to  both  hydrogen  electrode  and  oxygen 
electrode. 

At  the  interface  of  the  hydrogen  electrode  and  the  electrolyte 
(x  =  Jh2  ,  the  thickness  of  the  hydrogen  electrode),  the  rate  of  H2 
transport  to  the  surface  of  the  hydrogen  electrode  is  equal  to  the 
H2  production  rate  under  steady  state  conditions.  This  boundary 
condition  can  be  written  in  terms  of  the  electrical  current  density 
/  as: 


J 

2 F 


Pu2Bg  9Fh2 
RT/jlu2  dx 


(8) 


The  Dirichlet  boundary  condition  can  be  applied  at  the  sur¬ 
face  of  the  hydrogen  electrode: 


Rh2 


l*=o 


(9) 


Solving  Eq.  (7)  with  boundary  conditions  expressed  in  Eqs. 
(8)  and  (9),  the  pressure  of  H2  at  the  interface  of  the  hydrogen 
electrode  and  the  electrolyte  can  be  obtained: 


77?7/xh2Jh2 

PBg 


(10) 


Substituting  Eq.  (10)  into  Eq.  (3),  the  concentration  overpo¬ 
tential  of  the  hydrogen  electrode  can  be  expressed  as: 


JSOSE  =RL\n 

»conc,rl2  2  f 


VUh2)2  +  (JRTm2dH2/FBgi 


(11) 


H2 


At  the  oxygen  electrode  side,  the  mass  transfer  phenomena 
become  more  complex  as  both  diffusion  driven  by  concentra¬ 
tion  gradient  and  permeation  driven  by  pressure  gradient  occur 
simultaneously.  The  dusty  gas  model  (DGM)  is  employed  to 
describe  the  mass  transfer  within  the  porous  oxygen  electrode. 
Considering  both  diffusion  and  permeation,  the  DGM  can  be 
written  as  [18]. 


Ni 


neff  ner 

UUk  j=i  j#  Uij 


+ 


V  yj^i-yiNj 

neff 


1 

RT 


Pd  Vi  ,  d  P 

P - h  y  i  — 

dx  7  dx 


1  + 


hi  \ 

mj 


(12) 


where  DfJ  is  the  effective  Knudsen  diffusion  coefficient  of 
species  i;  D ??  represents  the  effective  binary  diffusion  coef¬ 
ficient  of  species  i  and  /;  yt  represents  molar  fraction  of  species 
i\P  is  the  local  total  pressure  in  the  oxygen  electrode;  and  /xm 
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is  the  dynamic  viscosity  of  the  gas  mixture  (O2  and  H2O).  The 
procedure  for  determination  of  molecular  diffusion  and  Knud- 
sen  diffusion  coefficients  can  be  found  in  literature  [12].  The  gas 
mixture  viscosity  (/xm)  can  be  obtained  by  Wilke’s  method  [19]. 

The  above  DGM  can  be  solved  numerically.  The  DGM 
domain  is  discretinized  for  the  finite  difference  method  (FDM) 
and  an  iterative  approach  with  considerable  sub-relaxation  is 
employed  to  obtain  the  distribution  of  gas  composition  and  total 
pressure.  Computation  is  repeated  until  convergence  is  attained. 
Substituting  the  calculated  partial  pressures  of  O2  and  H2O  at 
the  interface  of  the  oxygen  electrode  and  the  electrolyte  into  Eq. 
(4),  the  concentration  overpotential  at  the  oxygen  electrode  side 
can  be  obtained. 

3.  Results  and  discussion 

3.1.  Model  evaluation  and  temperature  effect 

The  above  model  has  been  derived  to  determine  the  J-V  char¬ 
acteristics  of  an  RSOFC-H  performing  dual  functions  of  fuel 
cell  and  electrolyzer.  The  model  evaluation  for  SOFC  mode  was 
provided  in  a  previous  study.  In  this  study,  the  newly  developed 
model  for  the  SOSE  mode  will  be  evaluated  by  comparing  the 
modeling  results  with  experimental  data.  In  the  literature,  only 
Iwahara’s  group  [7,9,10]  reported  the  detailed  J-V  characteris¬ 
tics  of  hydrogen  production  by  proton-conducting  based  SOSE 
and  thus  their  experimental  data  were  used  for  model  evaluation. 
In  these  experiments,  the  J-V  characteristics  of  an  electrolyte- 
supported  SOSE  cell  were  measured.  The  electrolyte,  hydrogen 
electrode,  and  oxygen  electrode  were  made  of  SrCeo.9  Y0.1 03_a, 
Pt  and  Pt  +  Lao.4Sro.6Co03_a  with  thicknesses  of  about  500  jxm, 
50  p,m,  and  50  |xm,  respectively.  The  J-V  data  were  mea¬ 
sured  at  temperatures  of  1073  K  and  1173  K  and  a  pressure  of 
1  atm. 

Fig.  3(a)  shows  the  comparison  between  the  electrochemical 
modeling  results  and  the  experimental  data  by  Iwahara  et  al.  [  10] . 
The  values  of  the  input  parameters  for  the  modeling  analysis  are 
shown  in  Table  1.  The  SOSE  potential  decreases  significantly 
with  increasing  temperature  (from  1073  to  1173K).  This  is 
mainly  because  the  Ohmic  overpotential  decreases  significantly 
with  increasing  temperature  due  to  higher  proton  conductiv¬ 
ity.  It  can  also  be  seen  that  the  simulation  results  generally 
agree  with  the  experimental  data.  The  relatively  large  discrep¬ 
ancy  between  simulation  results  and  experimental  data  may  be 
caused  by  inaccurate  values  of  electrolyte  resistivity  used,  which 
were  measured  from  Iwahara  et  al.’s  curves  [7].  In  addition,  the 
exchange  current  densities  used  for  calculating  the  activation 
overpotential  differ  significantly  in  the  literature  [11,13,20,21], 
which  could  be  another  source  of  inaccuracy.  In  order  to  better 
fit  the  experimental  data,  the  key  parameters  have  been  changed 
and  an  additional  comparison  with  experimental  data  at  a  tem¬ 
perature  of  1173K  is  presented  in  Fig.  3(b).  It  can  be  seen 
that  an  electrolyte  resistivity  of  2.2  £2  m  yields  better  simulation 
results  (Fig.  3(b))  than  2.0  Q  m  (Fig.  3(a)).  In  addition,  when  the 
exchange  current  density  is  changed  from  4000  to  200  Am-2, 
the  simulation  results  fit  the  experimental  data  even  better.  The 
findings  imply  that  the  simulation  results  can  fit  experimental 


Fig.  3 .  Comparison  between  theoretical  simulation  results  and  experimental  data 
by  Iwahara  et  al.  [10]:  (a)  Simulation  parameters  from  Table  1  and  (b)  Different 
simulation  parameters  at  1 173  K. 

data  well  by  adjusting  the  key  parameters,  such  as  electrolyte 
resistivity  and  exchange  current  density  of  the  electrodes. 

In  order  to  better  understand  the  characteristics  of  the  SOSE 
cell  for  hydrogen  production,  the  individual  Ohmic,  activation 
and  concentration  overpotentials  were  studied  as  shown  in  Fig.  4. 
The  results  show  that  the  Ohmic  overpotential  dominates  the 
activation  and  concentration  overpotentials.  At  a  current  density 
of  5000Am-2,  the  Ohmic  overpotential  is  above  8  V  (Fig.  4a), 
while  both  the  activation  and  concentration  overpotentials  are 
about  two  to  four  orders  of  magnitude  less  than  the  Ohmic 
overpotential  (Fig.  4b  and  c).  This  is  because  of  (1)  low  ionic 
conductivity  of  the  proton-conducting  electrolyte  and  (2)  thick 
electrolyte  (i.e.  500  |xm)  used.  The  overall  effect  explains  why 
the  J-V  curve  approximately  follows  a  linear  relationship.  The 
extremely  high  Ohmic  overpotential  also  rules  out  the  feasibility 
of  electrolyte- supported  configuration  for  RSOFC-H. 
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Table  1 


Values  of  input  parameters  used  in  the  present  study 


Parameter 

Value 

Operating  temperature,  T  (K) 

For  model  evaluation 

1173  and  1073 

For  parametric  analyses 

1073 

Operating  pressure,  P  (atm) 

1.0 

Gas  composition  at  the  surface  of  oxygen 

0.9/0. 1 

electrode  (molar  ratio  of  H2O/O2) 

Electrode  porosity,  s 

0.4 

Electrode  pore  radius,  r  (pun) 

0.5 

Electrode  tortuosity,  £ 

5.0 

Electrolyte  resistivity,  tp  (£2  m) 

For  model  evaluation  (SrCeo.6Yo.i03_a) 

3.3  (1073  K) 
2.0(1173  K) 

For  parametric  analyses  (SrCe0.95Y0.05  03-0- 

1.0 

or  BaCeo.8Smo.202.9) 

Electrolyte  supported  RSOFC-H 

Thickness  of  electrolyte,  L  (fxm) 

500 

Thickness  of  hydrogen  electrode,  dn2  (|Jim) 

50 

Thickness  of  oxygen  electrode,  do2  (|xm) 

50 

Hydrogen  electrode- supported  RSOFC-H 

Thickness  of  electrolyte,  L  (fxm) 

50 

Thickness  of  hydrogen  electrode,  du2  (|Jim) 

500 

Thickness  of  oxygen  electrode,  do2  (|xm) 

50 

Oxygen  electrode-supported  RSOFC-H 

Thickness  of  electrolyte,  L  (|xm) 

50 

Thickness  of  hydrogen  electrode,  du2  (|xm) 

50 

Thickness  of  oxygen  electrode,  do2  (|xm) 

500 

3.2.  Concentration  overpotentials  ofRSOFC-H  in  the 
SOSE  mode 

In  practice,  a  cell  consists  of  a  thick  layer  of  electrode  (typ¬ 
ically  500  pim)  that  can  provide  mechanical  support,  while  the 
other  electrode  is  a  thin  layer  (typically  50  p,m)  to  minimize 
the  overpotential  loss.  The  characteristics  of  different  electrodes 
(thick  or  thin,  hydrogen  or  oxygen)  are  shown  in  Fig.  5a.  A 
thick  hydrogen  electrode  only  causes  a  negligible  increase  in  the 
concentration  overpotential  in  comparison  with  a  thin  hydrogen 
electrode.  At  the  hydrogen  electrode,  the  reactant  (H+)  transport¬ 
ing  in  the  solid  electrolyte  is  not  limited  by  the  porous  structure. 
Therefore,  there  is  no  limiting  current  for  the  hydrogen  elec¬ 
trode.  In  addition,  the  small  molecular  weight  of  the  product 
(H2)  facilitates  its  transport.  As  a  result,  the  resistance  of  the 
porous  hydrogen  electrode  to  the  gas  transport  is  very  low,  lead¬ 
ing  to  low  concentration  overpotential.  For  comparison,  a  thick 
oxygen  electrode  suffers  seriously  from  high  concentration  over¬ 
potential  and  low  limiting  current  density  (around  4000A  m-2). 
Unlike  the  hydrogen  electrode,  the  reactant  at  the  oxygen  elec¬ 
trode  (H2O)  is  transported  through  the  pores,  and,  thus,  its 
transport  depends  on  the  porous  structures  of  the  electrode. 
The  high  molecular  weights  of  H2O  and  O2  (product)  result 
in  low  effective  diffusion  coefficient  of  H2O,  leading  to  high 
concentration  overpotential.  As  a  result,  a  hydrogen  electrode¬ 
support  performs  better  than  an  oxygen  electrode- support 
in  terms  of  low  electrical  energy  requirement,  as  shown  in 
Fig.  5b. 
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Fig.  4.  Simulation  results  of  individual  Ohmic,  activation,  and  concentration 
overpotentials  of  the  electrolyte-supported  SOSE  cell  based  on  proton¬ 
conducting  electrolyte  used  by  Iwahara  et  al.  [10]:  (a)  Ohmic  overpotential, 
(b)  activation  overpotentials,  and  (c)  concentration  overpotentials. 

3.3.  Concentration  overpotentials  ofRSOFC-H  in  the 
SOFC  mode 

The  concentration  overpotentials  of  an  RSOFC-H  working  in 
the  SOFC  mode  are  shown  in  Fig.  6a,  using  the  electrochemical 
model  from  Ni  et  al.  [1 1] .  Similar  to  the  SOSE  mode,  the  concen¬ 
tration  overpotentials  of  both  thick  and  thin  hydrogen  electrodes 
are  negligibly  small.  This  is  because  the  light  hydrogen  gas  is 
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Fig.  6.  Theoretical  modeling  results  of  electrochemical  performance  of  RSOFC- 
H  in  the  SOFC  mode:  (a)  concentration  overpotentials  of  electrodes  with 
different  thicknesses  and  (b)  power  densities  of  the  cell  with  different  support 
configurations. 


transported  in  the  hydrogen  electrode  quickly.  Similarly,  there 
is  no  limiting  current  density  for  the  hydrogen  electrode  in  the 
SOFC  mode.  For  comparison,  the  concentration  overpotentials 
of  the  oxygen  electrode  increase  considerably  when  a  thick  elec¬ 
trode  layer  is  used.  It  is  also  noted  that  the  limiting  current 
density  of  the  thick  oxygen  electrode  (500  p,m)  in  the  SOFC 
mode  is  about  1200 Am-2,  which  is  much  smaller  than  that 
in  the  SOSE  mode  (about  4000 Am-2).  In  the  SOFC  mode, 
the  molecular  weight  of  the  reactant  at  the  oxygen  electrode 
(O2)  is  higher  than  that  in  the  SOSE  mode  (H2O),  resulting  in 
lower  rate  of  transport  and  a  lower  limiting  current  density.  Due 
to  the  difference  in  concentration  overpotentials,  the  hydrogen 
electrode- support  performs  better  than  the  oxygen  electrode¬ 
support  in  the  SOFC  mode,  as  shown  in  Fig.  6b. 

3.4.  Comparison  between  RSOFC-H  and  RSOFC-O  and 
the  prospect  of  RSOFC-H 

From  the  analyses  in  Sections  3.2  and  3.3,  it  can  be  seen  that 
the  selection  of  an  electrode  support  can  greatly  affect  the  over¬ 
all  performance  of  an  RSOFC.  More  importantly,  the  hydrogen 
electrode- support  is  favorable  for  RSOFC-H  working  in  both 
SOSE  and  SOFC  modes.  For  comparison,  an  RSOFC-O  has 
challenges  with  selecting  either  the  oxygen  or  hydrogen  elec¬ 


trode  as  the  cell  mechanical  support  since  one  working  mode 
always  suffers  from  high  concentration  overpotentials.  There¬ 
fore,  the  RSOFC-H  has  a  natural  advantage  in  comparison  with 
RSOFC-O,  from  the  cell  structural  design  point  of  view. 

It  can  also  be  seen  from  the  present  study  that  the  Ohmic 
overpotential  of  RSOFC-H  is  yet  another  important  source  of 
energy  loss  because  of  low  conductivity  of  the  presently  used 
proton  conductors.  Due  to  its  high  Ohmic  overpotential,  the 
proton-conducting  SOFC  exhibits  lower  maximum  power  den¬ 
sity  (below  0.4  W  cm-2)  than  the  oxygen  ion  conducting  SOFC 
(above  0.6  W  cm-2)  [11].  In  order  to  further  enhance  the  per¬ 
formance  of  RSOFC-H,  it  is  highly  desirable  to  develop  novel 
materials  with  high  proton  conductivity.  Alternatively,  fabrica¬ 
tion  of  thin  film  electrolyte  can  greatly  reduce  the  electrolyte 
Ohmic  overpotential.  In  literature,  the  oxygen  ion-conducting 
electrolyte  of  conventional  RSOFC-O  can  be  made  as  thin  as 
8  Jim  [22] .  It  is  expected  that  the  hydrogen  electrode- supported 
RSOFC-H  can  be  further  developed  to  achieve  high  energy  con¬ 
version  efficiency. 

4.  Conclusions 

An  electrochemical  model  has  been  developed  to  investigate 
the  J-V  characteristics  of  an  RSOFC-H  based  on  proton- 
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conducting  electrolyte,  which  can  function  as  an  SOSE  to 
produce  hydrogen  fuel  and  as  an  SOFC  to  generate  electricity 
using  the  hydrogen  fuel  produced.  It  is  found  that  the  hydro¬ 
gen  electrode- support  is  favorable  for  both  SOSE  and  SOFC 
modes.  For  comparison,  a  conventional  RSOFC-O  based  on 
oxygen  ion-conducting  electrolyte  will  cause  high  concentra¬ 
tion  overpotential  in  either  one  of  two  working  modes  no  matter 
which  electrode-support  (hydrogen  or  oxygen)  is  selected.  The 
natural  advantage  of  RSOFC-H  makes  it  a  promising  technol¬ 
ogy  to  support  renewable  energy  utilization.  The  present  study 
suggests  that  development  work  on  electrolyte  materials  with 
higher  proton  conductivity  and  enhanced  fabrication  of  thin  film 
electrolytes  can  improve  RSOFC  performance. 
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